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ABSTRACT 

Using the nearly full sky K s = 11 .75 2MASS Redshift Survey [2MRS] of ~ 45, 000 galaxies 
we reconstruct the underlying peculiar velocity field and constrain the cosmological bulk flow 
within ~ 100 h _1 Mpc. These results are obtained by maximizing the probability to estimate 
the absolute magnitude of a galaxy given its observed apparent magnitude and redshift. At a 
depth of w 60 h *Mpc we find a bulk flow v B = (90 ± 65 , -230 ± 65 , 50 ± 65 ) km s" 1 in agree- 
ment with the theoretical predictions of the ACDM model. The reconstructed peculiar velocity 
field v that maximizes the likelihood is characterized by the parameter p = 0.323 ± 0.08. Both 
results are in agreement with those obtained previously using the ~ 23,000 galaxies of the 
shallower K s = 1 1 .25 2MRS survey. 

In our analysis we find that the luminosity function of 2MRS galaxies is poorly fitted by 
the Schechter form and that luminosity evolves such that objects become fainter with increas- 
ing redshift according to L(z) = L(z = 0)(1 +z) 



h2.7±0.15 
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1 INTRODUCTION 

Peculiar velocities arising from the cosmological growth of density 
fluctuations affect the estimation of the distances and luminosities 
of extragalactic objects from their measured redshifts, a spurious 
effect commonly known as redshift space distortions. The effect is 
expected to be systematic because of the large scale coherence of 
the peculiar velocity field in the ACDM model. In this respect, red- 
shift space distortions provide a unique tool to validate the ACDM 
and gravitational instability scenarios, to probe the underlying the 
velocity field and to estimate the growth rate of cosmological sim- 
ulation f{z) = ^77, where z is the redshift, a = (1 + z) _1 is the 
expansion parameter and D is the linear growth factor. The growth 
rate / m ainly depends on the mass density parameter Q. and the 
redshift (Peebles! 1 1980h . The exact dependence is determined by 
the underlying theory of gravity. Therefore the estimation of f(z) 
constitutes a sensitive test to Einstein's General Relativity. 

The availability of large and relatively deep galaxy redshift 
surveys in recent years has triggered a strong interest on the appar- 
ent anisotropy in galaxy clustering induced by redshift distortions 
since they can be used to tigh ten constraint s over different cos- 
mological parameters dAmendola et alj|2005l) and that to provide 
a unique way to discrimin ate between a dark energy scenario an d 
a modified gravity theory dGuzzo et afll2008l : Izhang et al1l2008l) . 



Redshift space distortions in the galaxy distribution at different 
epochs are now regarded as one of the most effective ways to attack 
the dark energy problem and constitutes one of the main scientific 
goals of generat ion redshift surveys like BigBoss 1 Schle gel et al.l 
l201ll) or Euclid dLaureiis et ai1l201 ll) . 

Redshift distortions not only misplace galaxies. They also af- 
fect the estimate of the objects' luminosities. And yet, little atten- 
tion has been given to this second aspect of the same phenomenon. 
In fact, the idea of using systematic biases in the estimated galaxy 
luminosities to constrain peculiar moti ons is not new. It dates back 
to the work of [jammann et al. (1979) who correlated the magni- 
tudes of nearby galaxies with their redshifts to constrain the veloc- 
ity of the Virgo cluster relative to the Local Group. However, the 
method requires a large number of objects to be effective. For this 
reason several authors focus ed on average quanti ti es rather than 
single objects. For exam ple iBaleisis et alj Jl998t> : iBlake & Walll 
d2002h : lltoh et al.l d2010h (and references therein) have exploited 
the Compton-Getting effect and searched for dipole variation in the 
surface nu mber density of dis t ant g alaxies to estimate the bulk flow. 
Similarly, lAbate & Feldmanl d201lh looked for a dipolar modula- 
tion in the variation of a suitably defined average apparent magni- 
tude across the sky. Currently available all-sky redshift surveys, like 
2MRS, allow an estimate of galaxy luminosities for a large num- 
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ber of galaxies. In this case, peculiar motions can be inferred from 
system atic variations in th e estimated galaxy luminosities across 
the skv. lNusser et al.l d201 lh adopted this approach and showed that 
i t can be used to est imate the the bulk flow in the local universe. 
dNusser et al.l d2012h took a step further and showed that the same 
idea could be used to constrain the growth rate of density fluctua- 
tions in the local universe, i.e. at z ~ 0. 

In this paper we use the new redshift survey of nearly full-sky 
2MAS S Redshift Survey ( 2MRS) of ~ 45000 galaxi es with K s < 
11.75 dHuchra et al.l |20 1 1[) to supersede the work of iNusser et al.l 
and Nusser et al. 120121) based on the previous release with 
a brighter apparent magnitude cut K s = 1 1 .25. The aim is twofold. 
First of all, there is a considerable interest in large-scale flows 
with some controversial claims of anomalous bulk fl ows on various 
scales t hat would exceed ACDM prediction (see e.g. lFeldman et al.l 
j2010h : lKashlinskv et al.l d2012l ) and refer ence therein), that wer e 
not confirmed by su bsequent ana l yses dNusser & Davis 
Bilicki et alj d201 lh: INusser et alj JEoTTh; iTurnbull et al 



2011); 



2012); 



Osborne et all d201ll) : iModv & Haiianl d2012l) . It is an intriguing 
issue that certainly justifies a closer look. Our technique provides 
a fresh approach to this outstanding problem. Thanks to the im- 
proved dataset we should be able to detect significant departures 
from ACDM predictions within 100 h _1 Mpc. Second of all, our 
technique constrains the velocity field independently on distance 
indicators. As such, it is free of potential systematic errors aris- 
ing from a miscalibration of the distance indicators. From the spa- 
tial distribution and estimated luminosities of 2MRS galaxies we 
are able to model the linear velocity field and determine its only 
free parameter, (3 = /(O.) /b, where b is the linear bias parameter 
of the galaxy sample. Our aim is to im prove the accuracy of the 
estimate obtained by INusser et alj J2012h and constrain the funda- 
mental quantity /(£2) at z ~ 0. 

Our method heavily relies on the estimation of the galaxy lu- 
minosity function [LF]. Therefore we take special care in detect- 
ing, evaluating and correcting for systematic biases related to the 
measurement of the LF. For this purpose we use a suite of differ- 
ent methods, some of which completely new, to estimate the LF 
and the selection functions of the sample. Since velocities are es- 
timated at the redshifts of the objects, our results are probe to the 
so-called Kaiser rocket effect. Any method aimed at estimating the 
underlying mass density field from a spatial distribution of mass 
tracers in a redshift survey should take this correction into account. 
Therefore, we perform this correction and, in Appendix, we offer 
an analytic treatment of the effect. 

The structure of the paper is as follows. In Section [2] we de- 
scribe the theoretical tools used in this work: the maximum like- 
lihood method and the different LF estimators used in the analy- 
ses. In Section [3] we describe the real and simulated datasets. In 
Section [4] we introduce, implement and apply a novel technique 
to estimate the selection function of the catalog directly from the 
observed redshift distribution of the galaxies. The results are com- 
pared with those obtained in Section [5] in which we first compute 
the LF and by integration, we obtain the selection function. In Sec- 
tions|6]and[7]we apply our maximum likelihood method to estimate 
the bulk flow and the p parameter. Finally, in Section[8]we discuss 
and summarize our main results. 



2 THEORETICAL TOOLS 

The scope of this work is to use the 2MRS flux limited at K s = 
1 1 .75 to estimate the cosmological bulk flow in the local uni- 



verse and trace the underlying peculiar velocity field. Peculiar ve- 
locities contribute to the measured redshift of an object. If r is 
the proper distance of a galaxy, z its redshift and v is the line 
of sight component of its peculiar velocity, then s = cz = r + v, 
where all quantities are expressed in kms , including the speed 
of light c. The absolute magnitude is estimated from the appar- 
ent magnitude, m, through Mq = m — 1 5 — 5 log cz. The "observed" 
magnitude is different from the true one M = m— 15 — 51ogr; = 
Mq — 51og(l —v/cz), where r; = r(l + z) is the luminosity dis- 
tance. The difference between M and Mo can be used to infer the 
peculiar velocity of the object. This can be done by maximizing the 
probability P(Mq\cz,V) of a galaxy having an observed magnitude 
Mq given its redshift cz and peculiar velocity v: 



P(M \cz,v) 



*(M) 



f-<L®(M)dM 



(1) 



where <1>(M) is the luminosity function, M; = W; — 25 — 5 log r/, the 
apparent magnitude limit of the catalog is m; = K„ = 1 1 .75 and the 
expression is valid as long as errors in the measured redshifts are 
small (a cz /cz < 1). 



2.1 Bulk flow from magnitudes and redshifts 

INusser et alj d201 ll) presented a simple method to measure cos- 
mological bulk flows by minimizing systematic variations in the 
galaxy magnitudes estimated from the observed redshifts. 

The method can be illustrated by the following example. Let 
us assume that the peculiar velocity field is characterized by a bulk 
flow vb and galaxies have a Schechter LF dSchech ter 1980): 



*(L) =0.41n(10)** 



l+a 



exp 



L 



(2) 



In the large distance and small redshift error approximation the 
probability in Eq.[TJis 

0.41n(10) ('|>y +a e- Z "'/ i * 
P(L |c Z ;v B ) = FTTZTT? i, \ ' (3) 

where Lq = (1 — 2\-^/cz)Lq and L/ = (1 — 2vg/cz)L/ and absolute 
magnitudes are related to luminosities through M = — 2.51ogL + 
const. 

The presence of a bulk flow systematically increases or 
decreases the estimated luminosity of a galaxy Lq. Therefore, an 
estimate of the bulk flow can be obtained by maximizing the prob- 
ability P(Lq\cz\ vb) with respect to vg. It is trivial to generalize this 
approach to a generic form of the luminosity function. 

This method can be seen as a generaliza t ion of the maximum 
likelihood approach proposed by lltoh et all d201fj|) . In that case 
the bulk flow was estimated from the apparent dipole anisotropy 
modulation in the s urface number density of galaxi es in the the 
SDSS-DR6 catalog, dAdelman-McCarthv et alll2008l) . To estimate 
the bulk flow that method requires angular positions, apparent mag- 
nitudes and the photometric redshift of the galaxies. The method 
proposed here requires more information (spectroscopic redshifts) 
but allows an estimate of the bulk flow from a differential quan- 
tity, the LF, rather than an integral one (the number density of ob- 
jects). As a result, th is method is mor e sensitive to bulk flows than 
the one proposed by lltoh et alj |2010). A similar method that does 
not use spectroscopic reds hifts but only appa r ent m agnitudes has 
been recently proposed by lAbate & Feldmanl d201 ll) . In that case 
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the authors looked for systematic, dipole-like vari ations in the ap- 
paren t magnitude of the LRGs in the SDSS survey feisenstein et al.l 
[2003)). 



2.2 (3 from magnitudes and redshifts 

iNusser et al] d2012h have extended the maximum likelihood tech- 
nique discussed above to constrain the full linear velocity field v. 
The method requires redshifts, angular positions and apparent mag- 
nitudes of galaxies. Galaxy positions are given in redshift space and 
are used to compute the linear velocity field as a function of (3. De- 
tai ls of the computation o f the linear velocity field can be found 
in INusser & Davisl dl994l) . Predicted velocities are used to com- 
pute distances and estimate the true absolute magnitude of the ob- 
jects. The best fit value of (3 is found by maximizing the probability 
P(Mq\cz,v($)) over all galaxies in the sample. 

We stress that the linear velocity field is predicted from the 
galaxy distribution in redshift space. Such a procedure is prone to 
the so-called Kaiser rocket effect, a systematic error induced by 
estimating the selection function of galaxies using redshifts rather 
than distances. In this work we explicitly correct for this bias. A 
detailed treatment of the effect can be found in Appendix A. 



2.3 Estimators of the galaxy luminosity function 

The measurement of the LF represents a key step in the maxi- 
mum likelihood methods outlined above. To guarantee an accurate 
measurement and to minimize possible systematic errors we have 
used different estimators for the luminosity function, <J>(M), that 
we briefly describe below. 

• l/Viviax estimator [<jy ]. This s imple non-parametric estima- 
tor originally proposed by Ischmidtl (1968) weights each object by 
the maximum observable comoving volume in which it can be de- 
tected. It is the only estimator among those we have considered 
that is sensitive to large scale inhomogeneities in the galaxy distri- 
bution. 

• STY estimator [<fts r ft]. Thi s estimator has been originally pro- 



posed by Sandage et al 



It assumes a Schecther form for 
4>(M) (Schechter 1980) and computes the best fit parameters by 
maximizing the product of the probabilities for galaxies to have a 
magnitude M given the observed redshift z (Eq.QJ: 

P, = IkP{Mi\zi). (4) 

For a Schechter form 

*(M) = 0.41n(10)**10°' 4 ( a+1 >( M *- M > 

x exp(-10 a4 ( M -- M >) (5) 

The free parameters are M* and a and the normalization <I>*. The 
latter does not concern us here. 

• Stepwise estimator [ <fta P ]. This method, originally proposed 
by Efstathi ou et al.l dl988l) . is non-parametric and has been derived 
from a maximum likelihood approach. The unknown luminosity 
function <t(M) is discretized into Nf, — 1 magnitude bins over the 
range M\ < ■ • • < Mtq b so that 

*(M) = <t>,- for M i+1 ^M>Mi . (6) 

Given this stepwise form for <t>(M), the probability in EqQ]becomes 



where we take M, + i — M; = AM = const and assume that the ac- 
tual and limit magnitudes, M and M;(z), fall into the bins i and j, 
respectively. The Nf, free parameters <!>, are estimated by maximiz- 
ing the product of single galaxies' probabilities. 

• Spline based estimator [^jp/]. The maximum likelihood 
method used to constrain the velocity field v requires a smooth LF 
in input. For this reason we introduce a new estimator in which the 
unknown LF is approximated by a smooth piecewise function 



4>(M) = qj{M) for M i+1 > M > M, 



(8) 



where q, is a third degree polynomial satisfying the boundary con- 
ditions qi(Mi) = <S, and g,(M, + i) = <I>, + 1 and defined such that 
its second derivative d 2 <&/dM 2 is continuous over the magnitude 
range M\ —Mp/ b . The coefficients of the splines can b e efficiently 
compu ted using the standard techniques described in IPress et all 
dl992|) . Splines can be integrated to estimate the denominator in 
Eq. [T]and the coefficients <!>,■ are determined by maximizing the 
product of probabilities P s = TljPj(Mj\zj) extended to all galaxies 
in the sample. Splines obtained from this procedure may produce 
a noisy LF especially at the faint and bight end due to the lim- 
ited number of objects. In order to suppress these spurious wiggles 
we maximize a function which is the sum of log[/> s ] and a penalty 
function which acquires very large negative values when the third 
derivative of the splines is large. This procedure efficiently sup- 
presses the wiggles and still yields a set of best fit coefficients <!>,■ 
based on maximum likelihood considerations. 



3 DATASETS 

We use the recently compiled 2MRS catalog that contains all galax- 
ies brighter than K s = 11.75 with measured spectroscopic redshift 
select ed from the 2MASS XSC catalog of nearly one million ob- 
jects dHuchra et al.l l20lH) . The catalog is 97.6 % complete and 
mostly unaffected by interstellar extinction and stellar confusion 
over the region \b\ > 5° for 30° < I 330° and \b\ > 8° otherwise. 
The total sky coverage is over 91 %. The catalog contains about 
43,000 galaxies and therefore represents a significant improvement 
over the K s ^ 1 1.2 5 red shift catalog of 23,0 00 objects used in the 
INusser et alj 1 201 H) and lNusser et alj (12012!) analyses. 



P{Mi\zi) = 



(Mi-Mj)<i>j+AMZ, k< j<i> k 



(7) 



Although the catalog extends out to s ~ 30,000 kms we 
restrict our analysis to objects at smaller distances. To minimize 
incompleteness for nearby objects we consider a semi-volume lim- 
ited sample that contains all galaxies with s > s C ut and galaxies 
with s ^ Scut that, if placed at s cut , would be brighter than the mag- 
nitude limit. We set s cut = 3,000 kms -1 . The semi volume lim- 
ited catalog is therefore obtained by excluding all galaxies with 
M > M[(s cut ) = m — 51ogr; (s cu t ) — 25 and s < s cut . In addition, 
throughout the paper we work with a version of the survey with 
collapsed fingers-of-god in the main nearby clusters and with the 
masked region near the galactic plane filled at a given redshift by 
folding the the positions of galaxies 

In addition to the real catalog we will also consider a suite of 
mock 2MRS catalogs. They will be used to test the validity of our 
likelihood approach and assess its uncertainties. Indeed we have 
used two different sets of mock catalogs: 

• Mock catalogs used to test the bulk flow accuracy 
[vb— Mocks]: 

These mock catalogs are the same as INusser et alj d201lh . The 
set is composed by 200 semi-volume limited mock catalogs with 
Scut = 3,000 kms -1 which contain the same number of objects as 
the real sample. Mock galaxies are randomly distributed within a 
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sphere of 200 h~'Mpc. Their absolute magnitudes are assigned ac- 
cording to the LF of 2MRS (early + late type) galaxies estimated by 
IWestovej {2007). Redshifts of the objects are the sum of the Hub- 
ble flow and peculiar velocities modeled as a random component 
sampled from a Gaussian distribution with zero mean and width of 
300 km s -1 . This scatter accounts for the combined effect of small 
scale velocity dispersion and errors in the measured redshift. No 
underlying large scale bulk flow was assigned to the mock galax- 
ies. 

• Mock catalogs used to test the ft accuracy [ft— Mocks] : 

These mocks are the same ones used in iNusser et al.l J2012h . 
They are a set of 135 2MRS mock catalogs extracted from the mock 
Two Micron All Sky Survey extracted from the Millennium simu- 
lation. Mock galaxie s were obtained the s emi a nalytic model of 
ISpringeletall j2005h : lDe~Lucia & Blaizotl J2007I) . In these mocks 
the central observer is not chosen at random. Instead, it is selected 
to match the density and the dynamical properties of our Local 
Group of galaxies. For our purposes the main relevant properties 
of the mock galaxies is their LF that is well approximated by a 
Schecht er form. More deta ils on the mock 2MRS galaxies can be 
found in lDavis et all fcOllh . 



4 REDSHIFT DISTRIBUTION OF 2MRS GALAXIES 
AND THEIR EVOLUTION 

The blue histogram in figure Q] shows the redshift distribution of 
objects in the semi-volume limited catalog of 2MRS galaxies, 
dN/ds. We use constant redshift bin of size As = 150 kms . 
Errorbars represent the Poisson scatter in each bin Gn = \/N = 
\J dN /ds x As. The black, continuous curve is a parametric fit to 
the distribution 



dN_ 
~ds 



■-AS{s)s 2 



(9) 



where S(s) is the galaxy selection function, i.e. the fraction of 
galaxies in the catalog at redshift s. The selection f unction can ei- 
ther be measured directly fro m the observed counts dKirshner et al.l 
ll979l : lDavis & Huchra|[l982T) or estimated from the LF. In this pa- 
per we adopt both approaches. In this section we use the first one 
and compare the result with the alternative approach in the next 
Section. For this purpose we present a novel method, the F/T es- 
timator, in which the selection function is computed by integrating 
the following equation 



dlnS(j) 



As- 



T(s) 



(10) 



where T(s) is the number of galaxies with redshift smaller than 
s that could also be detected at larger distances while F(s) is the 
number of galaxies within s that can only be detected out to s + As. 
The normalization A is set by matching the total number of galaxies 
in the catalog within s = 12,000 kms -1 . This estimator assumes 
that the galaxy luminosity and selection functions do not depend 
on the environment. 

In Fig. Q] we compare the smooth redshift distribution of 
galaxies obtained from the selection function estimated with the 
F/T method (black continuous curve) is compared with the ob- 
served dN/ds (blue histogram). The curve fits the data well out 
to s = 10,000 kms -1 . At higher redshifts the F/T method over- 
estimates the observed number of galaxies. If real, this difference 
would indicate that we are located within an under-dense region 
extending out to ~ 20, 000 km s~ 1 . Alternatively, the mismatch can 



observed 




s/H [Mpc/h] 



Figure 1. The blue histogram represents the observed redshift distribution 
of 2MRS galaxies computed in redshift bins A.v = 150 kms -1 . Errorbars 
represent 1 o~ Poisson uncertainties. Continuos lines show the fits to the 
data obtained from the selectio n function computed u sing the F/T estima- 
tor. The black curve assumes iKochanek et~aT] <200ll) k-correction but no 
luminosity evolution. The red curve shows the effect of k-correction a lu- 
minosity evolution e(z) = 3.04z. 



hint some systematic errors is affecting the data. Before consider- 
ing the first option, let us consider the possible sources of system- 
atic uncertainties. To estimate the ratio F(s)/T(s) in Eg. 1 101 one 
needs to compute the rest-frame absolute magnitude of galaxies, 
i.e to apply the so-called ^-correction, k(z). In addition, the galaxy 
luminosity may evolve with time according to some law, e(z). We 
correct for these systematic effects in our magnitude estimation as 
follows: 



M - 



■25-5logn{z)-k(z) 



(11) 



In the K s band and at low re dshift (z < 0.25) the k- correction is neg- 
ative. 0For 2MRS galaxies Kochanek et al. (2001), have found that 
k(z) = — 61og(l +z), a corrections that we adopted in this work. 
The evolution correction is more uncertain. It can be estimated by 
forcing a good match between the black curve and the histogram 
in Fig.[T] In practice we assume simple luminosity evolution model 
L(z) = L(z = 0)(1 +z) e and find £ by minimizing the y} function 



r 



Nbins 

: I - 



d7 



■As(sy 



(12) 



where the summation is over all redshift bins and a n is the Pois- 
son noise in the galaxy counts. We obtain L(z) = L(z = 0)(1 + 
z ~j+2.7±0.l5 jjj e resu i t Q f m j s correction is represented by the solid, 
red line in Fig.Q] which indeed provides a good fit to observations. 



5 THE LUMINOSITY FUNCTION OF 2MRS GALAXIES 

In this section we use the different estimators described in Sec- 
tion [5] to compute the LF of the 2MRS galaxies and compare the 



1 Since the K s band is on the Rayleigh Jeans part of the spectrum, and 
assuming the galaxy has little warm dust, k(z) is negative, unlike the optical 
bands. 
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Figure 2. LFs of 2MRS galaxies obtained from the different estimators de- 
scribed in the text. Q>v - blue triangles. Black, continuous curve. $>stp'- 
black dots. &spl'< rec '. continuous line. Errorbars represent 1-a bootstrap 
errors. Luminosity functions are computed in bins 0.25 mag. Symbols are 
plotted in every other bin for clarity. 
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Figure 3. Redshift distribution of 2MRS galaxies. Blue histogram: same 
dN/ds as in Fig.fT] Curves: dN/ds predicted from different LF estimators. 
Black, continuous : &sch- Red, continuos: &spl- Black, dashed 3>y. 



results to assess their robustness. Absolute magnitudes were com- 
puted in redshift space using the k(z) and e(z) corrections described 
above. LFs are estimated in bins of 0.25 mag. The results are plot- 
ted in Fig. [2] All LFs are normalized to match the number of 
galaxies in the bin at M = —23.37 which is close to the value of 
M* = —25.52 — 5 log h obtained from the Schechter fit. 

In the plot the blue triangles represent the LF obtained us- 
ing the 1 /Vmax method which is sensitive to large scale inhomo- 
geneities. There is an excess of faint (Mx s > — 20 + 5 log h) ob- 
jects that we interpret as an overdensity in the number of galax- 
ies in the nearby (s < 2,200 kms -1 ) volume of the universe. The 
black, thin curve shows the Schechter LF with best fit parameters 
(a = -1.0,M* = -25.52 -Slogft). At the bright end this fit de- 
viates significantly from the others. This is not surprising and re- 
flects that fact that we are forcing a Schechter fit to the LF of a 
composite sample of early and late type galaxies when we know 
that th e two subsample s are well fitted by two different Schechter 
forms dWestovej 2007) with M» for early t ype objects i0.5 m agni- 
tude brighter than for the late-type objects iNuss er et alj201 lb . The 
results of the <&stp (black dots) and <&spl (continuous, red curve) 
estimators are consistent with each other and with those of $>y for 
relatively bright objects Mk s — 5 log h < —20, i.e. in a volume large 
enough for fluctuations in the number density of galaxies to be neg- 
ligible. 

Errorbars in the plot are estimated from a bootstrap resampling 
analysis based on 20 random catalogs obtained from the original 
one by replacing each galaxy with n r objects, where n r is drawn 
from a Poisson distribution function with mean unity. All LF es- 
timators have been applied to the 20 bootstrap catalogs. The error 
bars represent the rms scatter among the catalogs. Bootstrap uncer- 
tainties turned out to be very close to Poisson errors. 

LFs are estimated in redshift space. We do not expect that they 
are different from those measured in real space, i.e. placing each 
objects at their true distances. The reason is that correction for pe- 
culiar velocities are of the order of ((v/cz) 2 ), where the average is 
taken over all directions in the surveyed area. Since we are dealing 
with an almost all-sky survey the effect is expected to be negligi- 



ble. To verify this hypothesis we have placed galaxies at their es- 
timated true positions using the peculiar velocity model discussed 
in Section [7] We found no significant difference between the LFs 
estimated in real and redshift space. 

From the estimated LF it is straightforward to compute the 
galaxy selection function: 



S{s) 



f M j {s) <J?(M)dM 



(13) 



for s > s cut . In Fig. [3] we show the dN/ds computed from the 
estimated selection function using Eq. [9] The different curves 
refer to different LF estimators. All of them are normalized to 
the observed galaxy counts (represented by the histogram) within 
12,000 km s" 1 . 

The dN/ds curves obtained from &sch (black continuous), 
$>Ste (not plotted) and $>s p i (red continuous) are independent of 
the underlying inhomogeneities in the galaxy distribution unlike 
the one obtained from <Ity (black, short-dashed) This difference is 
explains why the dN/ds obtained from the 1 /Vuax estimator is dif- 
ferent from the other curves both in the nearby region and in corre- 
spondence of the peaks of the distribution. The galaxy counts pre- 
dicted from the Schecther LF are systematically larger above the 
observed one at s > 10,000km s -1 . This discrepancy reflects the 
fact that the Schechter fit underestimates the number of bright and 
faint objects alike (Fig.[2j, resulting in an artificially small denom- 
inator in Eq . 1 1 3 1 

We notice the remarkable similarity between the dN/ds 
curves obtained from the spline and 1 /Vmcix estimators and the one 
obtained from the F/T method shown in figureQ] A result that fur- 
ther demonstrates the goodness of the new F/T estimator proposed 
in this paper. 

From the estimated LFs it is straightforward to compute the 
mean galaxy overdensity within s: 



1 + <§(<*)) 



j M : ( ' s) *(M)dM 



M,(s No , 



' ®(M)dM 



(14) 



The mean is set by the number density of galaxies within the red- 
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Figure 4. Mean galaxy overdensity (Y-axis) within the redshift s (X-axis) 
estimated from Q>sdi (blue, dashed curves) and &Spl ( re d> solid). Thick and 
thin curves indicate the two redshift cuts of 3,000 and 5,000 km s applied 
to extracted the semi-volume limited sample. All curves are normalized to 
(8(j< 12,000 kms" 1 )) =0. 



shift snorm = 12,000 kms chosen to normalize the selection 
function. The estimated cumulative overdensity obtained from the 
different LF estimators is plotted in figure [4] For the sake of clar- 
ity we only show the overdensity computed from the Schechter LF 
(thick, blue dashed line) and from the spline estimator (think red) 
since the latter is very similar to those obtained from the 1 /Vmox 
and the stepwise methods. The overdensity predicted by the two 
estimators are very similar within s = 12,000 kms~'. At larger 
distances, where the Schechter fit underestimates the number of 
expected counts, the corresponding overdensity is systematically 
smaller than obtained with the other estimators. These results de- 
pend on the choice of s cu , . Therefore we have repeated the analysis 
using a m ore aggressive cut s cu , = 5 , 000 km s~ 1 . The results of are 
shown by the thin curves in Fig. [4] Changing s cut only affects the 
estimate of the overdensity in the inner region 0< 3,000 kms" 1 ), 
as expected, but has no impact on the results obtained at larger radii. 



6 MEASURING THE BULK FLOW FROM LF 



We now apply the method described in Section [2T| to the 2MRS 
K s < 1.75 galaxy catalog to compute the bulk flow. Unlike 
iNusser et aLM20T l|) we do not break down the sample into early 
and late type objects since our best LF estimators ^t>spl an£ i &ste 
are designed to deal with a mixed population of objects. We take 
the LF estimated from ^spl as the reference case and will test the 
sensitivity of the results to the choice of the estimator. We do not 
considered the Schechter LF since, as we have shown, it underesti- 
mates the abundance of bright galaxies. 

We estimate the bulk flow in spherical shells of depth As = 
4,000 kms - ' and out to s = 10,000 kms - ' by maximizing the 
probability function P s = IT,P,(M()|cz, vb), where the product is 
over all galaxies and Pj is the probability of the single object de- 
fined in Eq. ??. T he binning in redshi ft and the volume sampled 
are the same as in INusser et al.l J201 lh but the number of objects 
is twice as large. The results are displayed in Fig. ??. Points con- 
nected by lines represent the Cartesian components of the differ- 
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Figure 5. The bulk flow of 2MRS galaxies with K s < 12.75 estimated 
from their luminosities in spherical shells 4,000 kms~'-thick. The three 
Cartesian components are represented by a dotted (X), solid (Y) and dot- 
dashed (Z) lines. 



ential bulk flow estimated at the mean radius of each redshift shell. 
Different line styles are used for the different Cartesian components 
specified by labels. Errorbars show the rms scatter in the bulk flows 
estimated from the 200 vb Mocks. These errors are contributed by 
shot noise, uncertainties in the observed magnitude and redshifts. 
They do not include cosmic variance. We find that the error budget 
is dominated by shot noise rather than redshift uncertainties.. Er- 
rors in the measured magnitudes induce the same systematic errors 
in both the measured and the reference LFs and do not affect the 
bulk flow estimate. Not surprisingly, the results are remarkably ro- 
bust to the method used to estimate the luminosity function: <t>spl 
and <l>s le and <fry . Finally, results are also robust to the ^-correction 
and/or the evolution corrections adopted, as we have verified by 
switching off either corrections and yet obtaining very similar bulk 
flows. 

Our results can be compared directly with those of 
dNusseretalj201ll) and shown in Fig. 1 of their paper. The two bulk 
flows are fully consistent with each other at all redshifts and for all 
Cartesian components. For example, at R ~ 6, 000 km s~ 1 we find a 
bulk flow of v B = (90 ± 65 , -230 ± 65 , 50 ± 65) km s" 1 , fully con- 
sistent with the one of dNusser et alj201 ll) : vb= (100±90, -240± 
90, ± 90)km s~ 1 . The fact that errors scale with the square root of 
the objects confirms that errors are dominated by Poisson noise. 



7 MEASURING P FROM LF 

In this section we apply the method described in Section 12.21 to 
model the linear velocity field traced by the 2MRS galaxies and to 
estimate the parameter p. The method uses the gravity field com- 
puted from the spatial distribution of the galaxies in the survey to 
minimize the scatter between the observed and the expected lumi- 
nosity function. The result provides an estimate of p. The gravity 
field is computed from the spati al distribution of galaxie s in redshift 
space using the linear method of lNusser & DavisUl994|) . As a refer- 
ence luminosity function we use the one determined with the spline 
estimator assuming the luminosity evolution and ^-correction dis- 
cussed in Section|4] Since the velocity field is found by maximizing 
P s = UjPj(Mo\z,\($)) with respect to p we need to produce a suite 
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of model v elocity fields for differ ent values of p. This is done by 
running the Nus ser & Davisl d 19941 ) reconstruction procedure using 
different values of P in the range [0.1,1] in steps of 0.02. 
The method is prone to several systematic errors: 

• Since galaxy positions are given in redshift space, a direct 
estimate of the mass overdensity would systematically affect the 
predicted peculiar velocities. We correct for this so-called Kaiser 
rocket effect as described in Appendix A. 

• To apply linear theory we need to filter out nonlinear contribu- 
tions. For this purpose we first smooth the galaxy distribution with 
a Gaussian window of radius 400 km s and then apply a second 
Gaussian filter of radius R s to remove residual nonlinearities. We 
choose R s = 600 kms -1 as a reference case and evaluate the ro- 
bustness of the results to the choice of R s . 

• Li near theory recovers the flow pattern reas onably well up to 
8^ 1. {Nusser et alii 199 ll : [Branchini et al. 2002). For this reason 
iNusser et alj J2012h removed all galaxies in regions with overden- 
sity above 8 CMf = 1 and verified that the results did not change 
significantly when a more conservative threshold 8 

cm — 2 was 

adopted. In our analysis we do not use any threshold and test the 
robustness of the results to overdensity cuts. 

Given the model velocity fiel d v(|3) and the e s timate d galaxy 
luminosity function we apply the INusser & Davis! j 19941) method 
to all galaxies within 10,000 kms -1 and minimize — lnP s = 
— J2,Tn/ > (Mo/|cz,-; v(p)) with respect to p. We do not consider ob- 
jects beyond 10,000 kms -1 due to the rapid decline in the ob- 
served number density of galaxies, in Fig. [6] we show the quantity 
A% 2 = — 21nP s (P) + 21nP J (P m ,-„), where P, mn is the best fit value 
of p found at the minimum of the curve. The width of the curve at 
A% 2 = 1 provides an estimate of the l-o error. When all objects are 
considered (black, continuous curve) we find P = 0.323 ± . 035. t o 
be compared with P = 0.35 ±0.05 found by INusser etail J2012h . 
The errors quoted here are mainly contributed by shot noise, as 
confirmed by the scaling with the number of galaxies. However, 
other sources contribute to the error budget. We obtain more real- 
istic errors by repeating the analysis on the 135 P— Mock catalogs. 
These errors include contributions from cosmic variance, nonlin- 
ear effects as well as shot noise. We find that cosmic variance 
and shot noise are similar in size and that the best fit estimate is 
P = 0.323 ±0.083. 

To detect possible systematic errors and check the robustness 
of the results we have repeated the likelihood procedure changing 
the free parameters of the model. 

• The method relies on the measurement of the galaxy LF. We 
have checked that the choice of the estimator has a negligible effect 
on the P value. Thi s resul t is not surprising since, as shown explic- 
itly in lNusser et alj J2012h the form of the luminosity function only 
affects the weighting given to galaxies in a certain luminosity range 
but has no impact on the value of the best fit p. 

• The choice of galaxy evolution model adopted e(z) might also 
affect the results. However, the effect is expected to be small since 
we only consider objects within s = 10,000 kms -1 where the ef- 
fect of the evolution is small, as shown in Fig. [T] To verify this 
assumption we have repeated the analysis by switching off lumi- 
nosity evolution, e{z) = 0. The results, shown by by the red dotted 
curve in Fig. [6] confirm that the the effect is negligible 

• To assess the impact of nonlinearities associated to the peaks 
of the density field we have repeated the a nalysis excluding a ll 
galaxies above 8 CU , = 1 and 8 C1 « = 2, as in INusser et alj (2012). 
We find no significant change in the best fit P value. 




0.3 0.35 0.4 

P 

Figure 6. A% 2 as a function of p. Solid line refers to all galaxies within 
10,000 kms with magnitudes corrected for luminosity evolution. The 
best fit value is p = 0.323 ±0.035, where the quoted error is the width of the 
curve at A% 2 = 1 and account for the finite number of galaxies. Red dashed 
curve: same galaxy sample but no correction for luminosity evolution. (P = 
0.319 ±0.033). Blue dashed: subsample with galaxies at s > 4,000 km s _1 . 
(P = 0.317 ±0.05). 

• The smoothing radius R s may also affect the estimate of p. In- 
creasing the smoothing radius decreases the amplitude of the grav- 
ity field so that a larger value of P is required to fit the data. We 
have estimated the amplitude of the effect by increasing the Gaus- 
sian smoothing to R s = 1,000 kms -1 . The best fit value of p in- 
creases, as expected. However, the amplitude of this shift is ~ 3 
times smaller than the random error. 

• The catalog incompleteness for faint objects preferentially af- 
fects the innermost part of the sample. To minimize possible sys- 
tematics errors we have repeated the analysis excluding all objects 
within s = 4000 km s~ 1 . The result is illustrated by the blue, dashed 
curve in Fig. [6] The main effect is to increase shot-noise errors (by 
~ 40 %) with a negligible shift in the best fit p. 



8 DISCUSSION AND CONCLUSIONS 

In this work we have used the new 2MRS galaxy catalog of galaxies 
brighter than K s = 1 1.75 to estimate the peculiar velocity field and 
its bulk flow in the nearby (s ^ 10,000 kms -1 ) universe from the 
apparent brightening/dimming of galaxy luminosities. To do this 
we have used the sam e ma ximum likelihood te chniques proposed 
bv lNusser et al. I J201 1) and Nusser et alj J2012D and applied to the 
brighter 2MRS K s = 1 1 .25 sample, which contains about ~ 50% 
less objects. Since our technique heavily relies on the knowledge of 
the luminosity and selection functions of the galaxies, most of our 
efforts have been devoted in measuring accurately these quantities 
to keep possible systematic errors under control. 

The main results of our analysis can be summarized as fol- 
lows: 

• We have computed the LF of 2MRS galaxies down to Mk s — 
5 log/? = — 17 using four different estimators to searc h for possible 
biases in the measured LF (see e.g. Illbertetalj ( l2004l) ) and derived 
quantities. The LF obtained from the 1 /Vmox estimator shows an 
excess of objects fainter than M% s — 51og/i = —20 not seen with 
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other estimators. We interpret this feature as large scale overdensity 
in the distribution of nearby (cz < 2, 200 km s~ 1 ) galaxies since the 
1 /Vmox method is sensitive to large scale inhomogeneities. 

Parametric estimators in which the LF is modeled with a 
Schechter form may induce systematic errors. More specifically, 
we find that the best fit Schecther form to the observed LF systemat- 
ically underestimates the number density of bright (Mk s — 5 log h < 
25.5) objects. This bias arises from fitting a single Schechter lu- 
minosity funcjipn_to acomposite sample of early and late type 
objects dNusser et alj 1201 lh . We have introduced a novel, non- 
parametric estimator fo r the LF similar to the stepwise method of 
lEfstathiou et al.1 |l98S) with the advantage of producing a smooth 
LF to be fed into the maximum likelihood method used in this 
work. Both the new estimator and the stepwise method hint at an 
excess of faint objects, quantified by a steepening of the LF in ex- 
cess of the typical power-law slope a = —1.0. The significance of 
this faint end excess, however, is barely above 1 o. 

It turns out that the differences among the LFs estimated with 
various method are rather small and have no practical impact for 
our goals. That is to say that the bulk flow and the (3 values ob- 
tained from the maximum likelihood procedure do not significantly 
depend on LF estimator adopted. 

• From the LFs we estimate the selection function of the cat- 
alog and the redshift distribution of the 2MRS galaxies. The re- 
sults are compared with the observed dN/ds and with the redshift 
distribution obtained from a novel F/T estimator. This new statis- 
tical tools allows one to measure the selection function of the cat- 
alog from the observed redshift distribution of the galaxies, rather 
than from the LF. The F/T estimator systematically over predicts 
the number of objects in the outer part of the survey that we in- 
terpret as an evolution in the galaxy luminosity. According to this 
interpretation the galaxy luminosity evolves so that objects would 
grow fainter with the redshift. A simple luminosity evolution model 
L(z) = L(z = 0)(1 +z) +2 ' 7±015 is sufficient to fit to the observed 
dN/ds. 

This result is not surprising since galaxy evolution in the infrared 
band can be significant even in the local patch of the universe. In 
the far infra red, evidence for an evolution in the number density 
of IRAS galaxies has been reported by a number of authors, al- 
though there is some controversy on the magn it ude of the effect ( 



(1998); 



Oliver et alj dl992h:lFisher et al.1 dl995l) ; ISpringel & White] 
Takeu chi et alJ d2003l)~ and reference therein). However, all 



these works find that the galaxy number density increase with red- 
shift whereas in our case the positive evolution in luminosity would 
decrease the number density of objects selected above a given flux. 
This puzzling result, however, does not affect the outcome of our 
maximum likelihood methods since in our analysis we consider 
galaxies with s < 10,000 kms~', where the evolution is found to 
be negligible. 

• Our estimation of the bulk flow in the local (s < 
10,000 km s -1 ) universe is fully consistent with that of 



iNusser etal] d201ll) . For example, at R w 6,000 kms 1 we find 
v B = (90 ±65, -230 ±65, 50 ±65) kms" 1 to be compared with 
v B = (100±90,-240±90,0±90)kms-'. In both cases errors 
are dominated by shot noise and which explains the ~ 40% im- 
provement in the accuracy of the estimate. Cosmic variance is 
not included. This the most precise estimate of the bulk flow 
in the local universe obtained without using distant indicators. 
Among alternative bulk flow estimators, th e one used here i s ar- 
guably more preci s e than those proposed by lltoh et all d20ld) and 
lAbate & Feldmanl d201ll) since, in our case, we look for system- 
atic variations in a differential quantity, the LF, rather than using 



int egral quantities. Our method is also superior to the one proposed 
by lHaehnelt & Tegmarkl d 19961) that exploits the kinetic Sunyaev- 
Zer dovich CMB distortions along the line-of-sight to galaxy clus- 
ters dSunvaev & Zeldovichl 1 19801) . due to the limited number of 
available galaxy clusters within 10,000 kms -1 . Indeed, the like- 
lihood method proposed here is the idea tool to probe bulk flows at 
different locations in the universe from luminosity variations mea- 
sured in deep and wide n ext generations reds hift surveys (for a 
quantitative assessment see Nus ser et al.l d201lh ). 

Our estimate of the bulk flow ag rees well with t he re- 
cent estimates of the bulk flow of INusser & Davisl d201lh 
v B (R ~ 6,000 kms -1 ) = (120± 40,-250 ± 40,40 ±60)kms~ 1 . 
base d on the SFI++ Tully-Fisher catalog of distance indica - 
tors dMasters et ai]|2006l : Ispringob et al.ll2007l ; iDavis et alj|201 lh . 
We also agree with the bulk flow obtained from SNe la 
data dColin et alj 1201 ll ; iDai et al.1 l201ll) and most notably with 
lTurnbulletal.ld2012l) who used the high quality 'First amendment' 
dataset. Finally, our result also agrees with th e recent analysis of 
the 2MASS galaxy dipole lBilicki et alj d2oTTh . It is quite remark- 
able that consistent bulk flows have been obtained from different 
datasets and using different methods affected by different systemat- 
ics and that all of them are in agreement with ACDM expectations. 

• Using the gravity field computed from the spatial distribution 
of 2MRS galaxies we were able to obtain a linear model for the 
velocity field in the local universe by minimizing the scatter in the 
LF. From this procedure we obtain [3 = 0.32 3 ± 0.083, in agree- 
ment with the results of INusser et ail d2012l) ((3 = 0.323 ± 0.1). 
The increase in the accuracy is more modest than for the bulk 
flow since in this case the error budget is almost equally con- 
tributed by shot noise and cosmic variance, which is now accounted 
for. We did check that this result is robust to luminosity evolu- 
tion, to the choice of the LF estimator, to catalog incompleteness 
for faint objects and to the smoothing scale used to remove non- 
lineari ties. Our es t imatio n of |3 is also in good agreement with 
that of IDavis et all d201ll) obtained by comparing the gravity field 
of 2MRS Kf < 11.25 galaxies with the peculiar velocities in the 
SFI++ galaxy catal og ((3 = 0.325 ± 0.045) and also with that of 
iBilicki et al.1 d201lh ((3 = 0.38 ± 0.04) obtained from the 2MASS 
galaxy dipole. 

From (3 it is possible to constrain the growth rate of density 
fluctuations /(O.) = (3/fo if b can be determined independently. 
However measuring the galaxy bias from existing datasets is dif- 
ficult and no theory of galaxy evolution is currently able to pre- 
cisely constrain its value. Alternatively, one could estimate some 
parameter combination that does not depend on bias and yet pro- 
vides the possibility of efficiently discriminating among different 
cosmological scenarios. The combinations /(fl)ag = (3og ga [ pro- 
vides such possibility since both (3 and Og^ a i, the rms density con- 
trast in galaxy number, can be measured from redshift surveys and 
/(Q)c? 8 is as good as f(Q.) to test alternative cosmological sce- 
nario s dWhite et alj2009l ; |Percival & Whitel2009l : ISong & Percivall 
2009). The validity of this probe relies on the hypothesis that at 
8 h~'Mpc, the scale at which one measures CJg^/, galaxy bias 
has the same value as at the much larger scales in which one 
measures (3. Estimates of f(Q.)ag from the apparent anisotropy 
in galaxy clustering have been obtained in the redsh i ft range z ~ 
[0.2,0 81 dPercival et al.l|2004 ISamushia etail |2012| ; iBlake et al.1 
|2011UGuzzo etalj 120081) . Several authors have pointed out that 
an accurate estimate of the growth rate at z ~ 0, which can only 
be obtained from the peculiar velocity field in the local universe, 
would significantly increase the discriminatory power of this cos- 
mological probe. The most recent estimates have been obtained 
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bvlDavisetall d201lh (/(£2)o 8 = 0.31 ± 0.06 at z < 0.033) and 
iTurnbull et aflfcOlA (/(£l)o g = 0.4 ± 0.04 at z & . 02). In our 
case, t aking Og iga ; = 0.97 ±0.05 from lWestoveil < l2007n : lReid et al.l 
d2010|) , we find /(£2)o 8 =0.31 ± 0.09 for z < 0.033, in good agree- 
ment with the other results. 

Since in many models of modified gravity the growth rate can be 
parametrized as f(£l) = fl Y (Linder 2005) many recent works have 
focused on comparing the estimated growth index y with the canon- 
ical value 6/11 of a ACDM universe. The value of y from our esti- 
mate of /( Q)qg must be consistent with those obta ined from the es- 
timates of iDavis et all fcOllMTurnbull et al.l J2012h . The value ob- 
tained combining the two estimates (y = 0.616±0.052 Hudson and 
Tu rnbull 2012, ApJ submi tted) is consistent with the one obtained 
bv lNusser&DavisI boill) from the bulk flow of SFI++ galaxies 
(y= 0.495 ±0.096) and with those obtained from galaxy clustering 
on larger scales. All of them agree with a ACDM standard gravity 
model. 
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where n„AV is the actual number of observed galaxies within AV 
and we have assumed that AV is so small that r,- = r. 

Since we are given the redshift s = r + u rather than distances, 
we can only compute S at the redshift position Sj of a galaxy. This 
will introduce a bias known as Kaiser rocket effect since it would 
induce a spurious component in the gravitational attraction esti- 
mated in redshift space. Note also that AV will change by the trans- 
formation to redshift space, which will amount to the usual redshift 
distortions. As a result, the estimated density contrast in redshift 
space is in |A3l 

, <, n 1 n a 1 /„ 1 d\nS \ , _ /, 1 d\nS 

1±5=-— = T?eT7 l + ~TT- " = 1+5 ' l + -TT- 11 
n S(r) n S[s) \ s dins J \ s dms 

(A4) 

where ^ -^Ly are directly determined from the data. The extra term 
involving the logarithmic derivative of 5 is the Kaiser term and can 
be evaluated at s instead of r when one assumes linear theory. Any 
method aimed at estimating the underlying mass density field from 
a spatial distribution of mass tracers in a redshift survey should take 
this correction into account. 
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A2 Kaiser effect for a volume limited survey 

It is often assumed that Kaiser correction is negligible in a volume 
limited survey. We show here that this is not quite the case. Le 
us consider a volume limited survey of all galaxies that are bright 
enough to be seen out to a distance r max , i.e. all these galaxies have 
M < Mi v = mi — 51og(r mflv ) — 15. The true density contrast at r < 

f max IS 



1±8, = ™ 



(A5) 



APPENDIX A: THE KAISER ROCKET EFFECT 

In this appendix we work out the expression for the so-called Kaiser 
rocket effect for two different types of catalogs: a flux limited and 
a volume limited redshift survey. 



Al Kaiser effect for a flux limited survey 

Assume we have a flux limited redshift survey of galaxies with ap- 
parent magnitudes < m;. The selection function S(r) is proportional 
to the average density and, apart from a normalization factor, is de- 
fined as 



S(r) 



M,(r) 



<£{M)dM 



(Al) 



where M; = m/ — 151og(r/) — 15, r\ is the luminosity distance in 
km s _1 and <5(M) is the galaxy luminosity function. Here, for con- 
venience, we drop the subscript ; from the luminosity distance. The 
mean number density of galaxies within r„ is defined as 



3 ^ _J_ 



(A2) 



I assume that S(r„) is sufficiently large to avoid using the other 
estimator which has J3 in it (that estimator minimizes the variance 
at the expense of biasing n). The true density contrast in a cell of 
size AV at r is then 



l + 8 r 



1 



' ieAV 



1 



n 1 

n Sir) 



(A3) 



where h is defined as IN j (431^,^) where N is the total number of 
galaxies in the volume limited survey. However, since we observe 
redshifts and not distances r = s — 11, objects are selected according 
to their absolute magnitude estimated in redshift space My = m — 
51og(s) - 15 w -5log{r) - 2. 17(w/i) - 15 andM = m- 151og(r) - 
15 ~ Mq + 2. 17(u/s). Thus, for example, in a volume element AV 
where u is positive, the physical limit M[ v is larger than its redshift- 
space estimate Mo./,,,. Therefore, for redshift space data, the mean 
number density within s is 



M h +2.VJu/s 



<S>(M)dM . 



(A6) 



where the integration extends out toM/,, ±2. \l{u/s) instead of M/ v . 
To correct for this effect the density contrast should be defined as 



J M J: <S?{M)dM 



1±5 



S(r lr , 



n S(r m 



no 
il 



^ 1 t/lnS\ u 
s d\nr I r r max 



(Al) 
(A8) 
(A9) 



where we have used the definition of S in Eq. lAll This is the anal- 
ogous of the Kaiser correction for a flux limited survey except that 
11 is multiplied by a quantity that is defined at r max rather than at r. 
For the 2MRS K s < 11.75 catalog, dlnS/dlnr grows from 0.92 at 
= 3,000 kms to 1.42 at 5,000 kms~' and the corresponding 
Kaiser term in a sub-volume limited catalog cannot be neglected. 
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